Fractal design for efficient brittle plates under gentle pressure loading 
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We consider a plate made from an isotropic but brittle elastic material, which is used to span a 
rigid aperture, across which a small pressure difference is applied. The problem we address is to 
find the structure which uses the least amount of material without breaking. Under a simple set 
of physical approximations and for a certain region of the pressure-brittleness parameter space, we 
find that a fractal structure in which the plate consists of thicker spars supporting thinner spars in 
an hierarchical arrangement gives a design of high mechanical efficiency. 

PACS numbers: 46.25.Cc 



I. INTRODUCTION 

In the construction of buildings and implements, opti- 
mising mechanical efficiency (minimising the amount of 
material used) has been practiced since time immemo- 
rial. Striking examples may be seen in the flying but- 
tresses and remarkably thin vaulted roofs of medieval 
European cathedrals. The underlying principles for re- 
ducing the mass of masonry were developed much later 
however, starting with the concept of thrust lines 

After Euler studied the buckling of struts, it was 
recognised that compression members are in general con- 
siderably less efficient than tension members. This line 
of reasoning waspursued in the twentieth century by re- 
searchers (e.g. [H H|) who also considered the question 
of couplings required to attach different parts of engi- 
neering structures to one another. Because of the scaling 
properties of Euler buckling and the cost in material for 
couplings of tension members, it is generally advanta- 
geous to have few compression members and many long 
tension members. An efficient structure such as a space 
frame will therefore typically resemble a tent in this re- 
gard; a pattern which is also seen in the anatomy of land 
animals where compressional limb bones are sheathed in 
a finely divided web of tendons and musculature || . 

In more recent years, the main focus in structural op- 
timization has switched to computational approaches. 
Various classes of problem are considered: for example 
under a fixed set of loads and using a fixed total volume 
of material, to minimise the compliance (energy storage) 
or the maximum stress in the structure. For trusses, the 
"ground structure method" [|| is typically applied, where 
the initial structure is specified by a set of n points, and a 
spaceframe described by (a subset of) the complete graph 
K n which has these points as nodes. The cross sections 
of the beams are then varied (potentially down to zero 
size). 

For the optimization of solid components, the naive ap- 
proach of drilling holes to reduce weight has been brought 
to a high degree of refinement with methods such as the 
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"SIMP" schemeffor "Solid Isotropic Microstructure with 
Penalization" ) 0, H[ where voxels of the material are al- 
lowed to have grey-scale values during the optimization 
process. These so-called topology optimization meth- 
ods allow for any number and shape of holes to emerge 
from the simulation, except that a minimum length scale 
(larger than the voxel size) is imposed 0] because of the 
tendency for optimal structures to contain many fine ten- 
sion members [i| which are difficult to manufacture. 

As noted above, the tendency towards fine subdivision 
of tension members is well known. However a rather dif- 
ferent possibility has emerged from two widely different 
fields. This is the idea that fractals 0, which are well 
known to have interesting mechanical properties when 
they occur in colloidal floes [l(J [HI Q3 and have recently 
been shown to give highly efficient transport networks 
[l3| might also occur as optimum mechanical designs: 

Firstly it has been noted that trabecular bone has a 
fractal architecture, which has been hypothesized to be 
related to mechanical efficiency 14[ . Unfortunately, more 
recent work indicates that this complex microstructure 
may not be the result of such a simple optimization [l5j . 

Secondly, in paleontology the complex suture patterns 
of ammonites have been seen as an adaptation for greater 
strength [l(| [ijj ; but again more recent work suggests 
that fractal morphology is not correlated to life at greater 
depth [TU . Note however that this evidence is not conclu- 
sive, since strength and efficiency are different quantities; 
greater efficiency from fractal design may be demanded 
by lower availability of minerals and (as a very tentative 
extrapolation from results in this paper) be possible only 
at shallower depth. 

Given the tantalising possibilities and the complexity 
of the systems previously studied, it would be interesting 
to construct a problem in mechanical structure optimi- 
sation which is simple enough to analyse in depth, yet 
encompassing enough physics to be non-trivial. Careful 
analysis might then suggest whether fractal design princi- 
ples may indeed give highly efficient structural solutions. 

To this end, we consider the following problem: sup- 
pose we have an aperture in a rigid structure, which we 
wish to span with a thin plate of an elastically isotropic 
but brittle material (Fig.[T]). Furthermore, suppose that 
we want this opening to support a pressure difference, 
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Thin plate 




FIG. 1: Cutaway view of a thin pressure-supporting plate 
which spans a square aperture, and is clamped at the edges 
of this to a rigid support. 



so that the spanning plate has a tendency to bow or to 
break. We aim to find the minimum amount of material 
required to support this pressure without breaking. We 
also impose the condition that the plate is unstressed 
in the absence of the pressure loading. Such a situa- 
tion might arise in different contexts; for example in the 
manufacture of silicon membranes for membrane emulsi- 
fication [l9[ or the construction of pressure vessels. 



II. PRECISE STATEMENT OF THE PROBLEM 

Consider a thin plate of elastic material with Young 
modulus Y and Poisson ratio v. We assume that when 
the plate is not subject to a pressure load, it has no 
quenched-in stresses, it lies in the x — y plane, and has a 
possibly non-uniform thickness b(x,y). 

The plate is attached to the boundary of a rigid aper- 
ture, so that even under load, the out-of-plane deflection 
and its gradient are zero at the boundary. We refer to 
these constraints as "cantilever boundary conditions" . 

Let Lq be a typical unit of length (for example of the 
same order of magnitude as the size of the aperture), 
then we assume that the plate is thin, in the sense that 
b/Lo <C 1. Let us suppose that the plate is now subject to 
a small pressure difference P between its two sides, and 
we define a non-dimensional pressure loading parameter 
p = P/Y. 

In order to state the problem, we first need an expres- 
sion for the elastic energy of the plate under such pressure 
loading. We are furthermore interested in the case where 
the deflection of the plate may not be small compared 
to its thickness, and we will therefore have to consider 
stretching or shearing of the middle plane of the plate 

Under these circumstances, a useful approximation is 
that of von Karman [20|, l2llj . A simple exposition may 
also be found in Ref. [22j, from which we take (with a 
change of notation) the results we need below. 

Consider a region of the middle plane of the plate near 
a point which has co-ordinates before loading given by 
(xo^yojO). Under a load p, the plate will deform, after 
which we define the new co-ordinates of the point under 



consideration to be 

(x + u(x ,y()),yo + v(x Q ,y Q ),w(x ,y )) , 

where u, v and w are assumed to be small and we have 
allowed the possibility of out-of-plane deformations w. 

The two-dimensional strain tensor of the middle layer 
after deformation is given by e with components 
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For a uniformly thick plate (or a plate with sufficiently 
slowly varying thickness [23j ) the equilibrium deforma- 
tion of the middle plane can be obtained by minimising 
the e nerg y U of Eq. |(TJ) with suitable boundary condi- 
tions 122], where 



U = U S + U B - Y / dxdy(pw) 



(1) 



and the stretching and bending energies of the plate can 
be simply derived from the corresponding expressions in 
Ref. 223 as: 



Us = J dxdy ™ [u [Tr(e)] 2 + (1 - ^)Tr(e 2 )} 

(2) 

dxdy 2A(l-^) { [Tr(i/)]2 ~ 2(1 ~ ") det ( H )\ 



where H is the Hessian matrix 
H(x,y) = 



d w d w 
dx 2 dxdy 
d 2 w d 2 w 



(3) 



(4) 
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In addition to this, we note that the maximum tensile 
strain experienced by the material is given by 



q eq 



d 2 x w dl w 
dl y w d'w 



(5) 



where q T = (cos 8, sin 6). In this expression, the first 
term comes from stretching from in- and out-of-plane 
deformation of the middle plane and the second from 
stretching produced by bending of the plate (which is 
maximal at a distance ±6/2 from the middle plane). 

In what follows, we assume that the material is brittle, 
in that it will fail if it experiences a tensile strain that 
exceeds a small number e <C 1. We therefore have the 
no-breaking condition: 



< e 



(6) 
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We can now state precisely the problem of finding a 
mechanically efficient plate: for a fixed aperture ge- 
ometry, and for each pair of values (e,p), we wish to 
find that function b(x, y) which minimises the integral 
/ b(x,y)dxdy, while at the same time achieving mechan- 
ical equilibrium (minimising the energy of Eq. ([1])) and 
subject to the condition of no breakage in Eq. (J6j) . In fact, 
this problem as stated is not well-posed; we therefore 
choose the extra restriction (discussed in section IIV A[) 
that |V6| < G everywhere, for some fixed number G. 

It is tempting to approach this problem by choosing 
an arbitrary function b(x,y), minimising the energy, and 
then locally removing or adding material according as 
whether e max is greater than or less than e. This is a 
possible strategy, but based on the calculations to follow 
we suspect that starting from a uniform b, this strategy 
will only produce a local minimum in the amount of ma- 
terial used. 



III. ONE DIMENSIONAL CASE 

A. Solution for cantilever boundary conditions 

Let us consider the case where we have a plate which is 
infinitely long in the y-direction. Furthermore, we impose 
the strong restriction that & is a function only of x, the 
direction across the plate. This will turn out to be very 
sub-optimal, but it is an instructive case for developing 
the argument. Under these conditions, v = 0, while u, w 
and b are functions only of x. The energy per unit length 
in the y-direction of the plate under pressure becomes 



U 
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(7) 



which must be minimised subject to the no-breakage con- 
dition 
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As stated above, the material is brittle in the sense 
e <C 1, and we want the plate to be thin (b/Lo <C 1). 
We also expect u/L$ -C 1, w/Lq <C 1, and from the 
geometry, 0(u) = 0(w 2 ), which we write for brevity as 
u ^ w 2 . 

Now, consider first the case where b is a constant and 
the plate is a strip given by x £ (—a, a). Let £ = x/a then 
the Euler-Lagrange equations obtained by minimising the 
energy of Eq. ([7]) are 
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(10) 



where p = (1 — v 2 )p and £ is the constant of integra- 
tion of the Euler-Lagrange equations which is still to 
be determined (physically, ( has the interpretation that 
Qb 2 1 (12a 2 ) is the strain experienced by the mid-plane of 
the plate). 

Noting that w = d x w = at x = ±a, and that w is 
symmetric about x = 0, we find 
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C sinh ( 



2Ccosh(«) 



(11) 



C sinh £ 
sinh(2C0 - 2CC 



4C sinh^ C 



(12) 



The condition which determines £ is that u — at x 
±a, which gives 



1 
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r;2 



2 coshC 
C sinh ( 
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sinh(2C) - 2C 



4C sinh^ C 



(13) 

The Euler-Lagrange equation, Eq. © for u shows that 
the strain produced by stretching at the middle plane 
is uniform over the plate, and therefore the no-breakage 
condition is 



2 h 2 



C 2 b 
12a 2 
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2a 2 



d 2 w 



< e. 



(14) 



1. Breakage dominated by stretching: £ large 

Consider first the case £ large, which corresponds to 
pa 4 /b 4 ^> 1, then bending occurs mainly in thin "bound- 
ary layers" near the edges of the plate (i.e. close to £ = 1 
and £ = — 1) which have a width A£ ~ l/£. In this limit, 
we find 



(288) 1 V /3 



(!) 



4/3 



(15) 



From the no-breakage condition of Eq. (fl"4"|) , we can cal- 
culate the required thickness of the plate. We find that 
bending and stretching both contribute significantly in 
the boundary layer (but bending is unimportant else- 
where) , and the necessary minimum thickness 6 m j n is 



2.885pe 



(16) 



It is also interesting to know £ in terms of p and e for the 



case when b = b D 



in order to estimate the thickness of 



the boundary layers. The result is 

C w 0.626e 2 p-\ 

Lastly, consistency with the fact that pa 4 /b 4 3> 1 implies 
that p < e 2 . 
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2. Breakage dominated by bending: ( small 

Next, let us consider the opposite extreme, where 
breakage is dominated by bending everywhere; i.e. £ is 
very small, equivalent to pa 4 /o 4 <C 1. This leads to 



c 



There are no boundary layers, and the expression for w 
becomes 

--f (!)*('-«■)*• 

The no-breakage condition is therefore: 
( 2 b 2 2pa 2 



12a 2 



b 2 



< e, 



where the first term corresponds to strain from stretching 
of the middle layer (both in-plane and out-of- plane) , and 
the second term is strain from bending. We find that 
breakage is governed only by the bending term in this 
limit, and 



2 l/2-l/2 £ -l/2_ 



Consistency with pa 4 /b 4 <§C 1 leads to p ^> e 2 , and thin- 
ness of the plate (b/a) <C 1 also imposes the condition 
p < e. 



B. Freely hinged boundary conditions 

In the case of £ large, the required minimum thickness 
of a uniform plate is determined by the thickness needed 
to sustain the bending and stretching forces in the thin 
boundary layers at the edges of the plate. Given more 
freedom to design a plate with non-uniform thickness, 
one might imagine that the most efficient design would 
have a larger value for b in the two boundary layers near 
the edges, and be thinner everywhere else. The required 
thickness of this middle region could be estimated con- 
servatively from the minimum thickness required for a 
plate with freely hinged boundary conditions. This case 
is also immediately solvable: we require that w be sym- 
metrical while u = 0, w = and d 2 w = at x = ±a. 
The Euler-Lagrange equations Eq. (J9]) and (fTT))) then give 
in the limit of large £ 



0.408pe 



-3/2 



(17) 



The important thing to note is that although the pre- 
factor is smaller in Eq. (fT7|) , the order of magnitude of 
all quantities is the same as for the cantilevered boundary 
condition case (Eq. (|16[) ). Even if we were to optimise 
the mechanical efficiency of the plate by choosing b to be 



thicker in the boundary layer and thinner in the centre of 
the plate, we would only achieve a gain in efficiency by a 
numerical factor of order unity. We also conjecture that 
an optimisation strategy comprising removing material 
where the breakage condition is not satisfied will converge 
on a solution of this kind. 

The boundary layers we encounter in this case are 
therefore in some sense unimportant; having only a quan- 
titative and not a qualitative effect on the amount of ma- 
terial needed to make our pressure-bearing plate. Shortly 
however, we shall encounter cases where the presence of 
boundary layers (regions where quantities vary rapidly 
in space) alters even the scaling of mechanical efficiency 
with p and e. 



C. Order of magnitude behaviour for a uniform 
plate 

We can now state the order-of-magnitude behaviour of 
various quantities (assuming v ^ — 1) for our optimum 
thin plate which has either uniform thickness, or bound- 
ary layers at the edges. In what follows, we shall use 
"~" to denote that two quantities are the same order of 
magnitude. We have: 



pc 



-3/2 



p<^c z 



pl/2 £ -l/2 £ 2 ^ p ^ £ 



(18) 



where the first case corresponds to a stretching- 
dominated and the second to a bending-dominated 
regime. 

The maximum deflection, maximum curvature from 
bending and maximum strain from stretching (both in- 
plane and out of plane) are given by 



w(0) 



,1/2 



p < e 



p- 1 / 2 e 3 / 2 e 2 < p < e 



d 2 w 
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-l 



P < € Z 

p < e 2 
e 2 < p < e 



in b. layer 
out b. layer 



2l2 



p < c z 



(19) 



(20) 



(21) 



These relations are shown in the top panel of Fig. [2] as an 
"order of magnitude" ("o.o.m.") plot. The advantage of 
the compressed scale used in this plot is that all crossover 
behaviour is sharp, and prefactors are irrelevant (indeed 
even a factor as large as log(e _1 ) would be invisible on 
this scale). The bottom panel of Fig.[5]shows the different 
regimes in the p — e plane. Because of the different scale, 
pre- factors are now important, and so the picture must 
be regarded as schematic (except very far from the origin 
at (loge.logp) = (0,0) ). 
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FIG. 3: Schematic for the design of a plate consisting of par- 
allel spars separated by long narrow panels, and spanning a 
square aperture of side 2L. 

Gs in Fig. [3]) the supporting spars are designed to be. 
However, at some point, very narrow and tall spars are 
likely to buckle under loading. In the analysis that fol- 
lows, we have therefore placed a simple restriction on the 
design of the plate, which is both interesting in its own 
right, but can also be seen as an attempt to capture this 
buckling limitation. The design restriction we place is 
that the maximum gradient of the thickness cannot ex- 
ceed some set value G, which is of order unity. There is 
also a third reason for choosing such a G, which repre- 
sents a limitation of the current analysis, rather than a 
change in the qualitative beahviour of the system. This 
limitation is that if the thickness b varies too quickly in 
space, then Eqs. @ and ([3]) will no longer be valid. 
Stated formally the new restriction is that 

V(x,y):\Vb(x,y)\<G. (22) 



FIG. 2: Top figure: order of magnitude behaviour for quan- 
tities as a function of pressure, for an aperture spanned by a 
plate of uniform thickness. Plotted on this scale, all crossovers 
are sharp and all order 1 prefactors are invisible. Bottom fig- 
ure: schematic picture of the p — e plane, showing the approx- 
imate positions of crossovers in behaviour for this case. The 
dashed line showing the crossover of regimes is p ~ e 2 



IV. PLATE WITH PARALLEL SPARS 
A. Behaviour of a single spar 

So far we have analysed the case where the thickness b 
is constant, or (by arguing from these solutions), the case 
where there are narrow boundary layers at the edges of 
the plate accommodating the local bending forces. How- 
ever, because the bending term (Eq. ([3J) in the energy 
(Eq. ([1])) has a pre- factor cubic in the plate thickness, 
we expect intuitively that it would be more efficient still 
to design a plate with narrow, thick (large b) beams or 
"spars" , which support thinner "panels" between. A sim- 
ple example of such a design would be as shown in Fig. [3j 

Again, from the form of the energy expression, we 
would expect that the design would get more efficient 
the narrower (smaller s in Fig. [3]) and the taller (larger 



Let us therefore consider a spar which has a diamond- 
shaped cross section satisfying this constraint, with width 
s, length 2L and loaded under a force / per unit length. 
We define a re-scaled parameter with dimensions of 
length to describe this loading, namely / = f/Y. 

In order to write down the energy of this system from 
Eq. ([1]), we need a further condition on the deformation 
of the middle plane. Because we are envisaging a narrow 
spar, we take this condition to be that the component 
of the stress in the x-direction (across the spar) at the 
middle plane is zero. Because the x and y directions are 
still the principal directions for both the strain and stress 



tensor, then this leads to e ra 
of the spar is given by 



-ve yy and so the energy 
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-(dyw) 2 



fw 



(23) 



Defining rj = y/L, the Euler-Lagrangc equations for Eq. 
(1231) are 



C 2 GV 



dv 

drj ~ 48L 2 (1 
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2L 
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G 2 s 2 d 4 w _ C 2 G 2 s 2 d 2 w 2/(1 - v 2 )L 2 
48L 2 drj 4 ~ A8L 2 !hf + 



Gs 2 



where £ is a constant of integration in the Euler-Lagrange 
equations for Eq. (f2"3")l and needs to be determined. Solv- 
ing these with cantilever boundary conditions gives 
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4608L 7 (1 



,2\2 f2 



y 



48L(1 - v 2 ) 1 
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( sinh £ 



G 6 s 8 C 4 

2 sinh(Cr/) sinh(2(>7) - 2(r] 



Q 2 sinh C, 



4Csinh 2 C 



where the constant of integration £ can be determined 
from 



G 



8 o 10 



2\3 



L S f 2 

1 2coshC 2 
3 ~ CsinhC + C 2 



221184(1 -v 1 ) 
sinh(2C) - 2( 



4C sinli C 



From now on, we consider only the orders of magnitudes 
of different quantities, and again there will be two cases: 



1. Stretching dominated regime 

In the stretching dominated regime, we have C 3* 1, 
and 



C' 



£4/3 jl/3 

G 4 / 3 s 5 / 3 : 



so / ^> G 4 s 5 j ' L A and the no-breaking condition (which is 
dominated by stre tching everywhere except at the ends 
of the spar, where bending is also significant) gi ves 

/ < e 3/2 s 2 L~ 1 G. 

Because of the no-breakage condition, the strain from in- 
plane and out-of-plane stretching is always less than e. 
The maximum amount of bending is 



d 2 w 



dy 2 



L 2 f L 2 / 3 f 2 / 3 



G 3 s 4 C G 5 /3 s 7/3 ' 



which occurs in the "boundary layers" near the two ends 
of the spar. The maximum out-of-plane deflection is 



w(0) 



£4/3 j 1/3 

GW/s- 



2. Bending dominated regime 

In the case that £ is small, we are in the bending dom- 
inated regime and we find that 



c 



G 4 s 5 



Therefore consistency requires that / -C G 4 s 5 /L 4 and 
the no-breaking condition (which is dominated by bend- 
ing) is 

/ < es 3 L- 2 G 2 . 
The maximum amount of bending is given by 



d 2 w 



dy 2 



L 2 f 
G 3 s 4 



and the amount of in-plane and out-of-plane stretching 
is of order 

( 2 G 2 s 2 L 6 f 2 



L 2 G 6 s 8 ' 
The maximum out-of plane deflection is 

B. Parallel spars supporting panels 

Let us suppose we have a plate consisting, as in Fig. [3l 
of parallel "spars" separated by thinner "panels" . The 
parallel spars also represent a collection of boundary lay- 
ers for the function b(x,y), in the sense that quantities 
have a fast variation as a function of a: or y. In contrast to 
the case discussed above, these boundary layers produce 
a qualitative change in the system's behaviour. 

We assume that the entire aperture is a square of side 
2L, and that the spars are of width s, and separated 
by a distance 2a from neighbouring spars. Because the 
analysis from now on will focus on orders of magnitude, 
we shall assume that G is of order unity, and omit it 
from the expressions. We will assume (and later check) 
that s <C a and a <C L (the latter, together with the 
relative rigidity of the spars, ensuring that we can treat 
the panels as infinitely long strips). The panels have a 
thickness bo- 

We now calculate the average thickness (b) of the plate 
that is required to support a pressure load, where (again 
in order-of-magnitude terms) 



(b) 



ab 



max(a6o, s 2 ) 
max(a, s) 



(24) 



and the second equivalence follows from the fact that 
the terms being added (in each of the numerator and 
denominator) will in general be of different magnitude. 
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TABLE I: The six conditions which must be satisfied by spars 
and panels for a plate consisting of parallel spars supporting 
intervening panels. Spars and panels may be either in the 
bending or stretching dominated regimes. 

Spars in stretching Spars in bending 

dominated regime dominated regime 

b ~ ape _3/2 b ~ ape~ 3/2 

Panels in p < e 2 p < e 2 

stretching (f ) « (£) 1/5 p 1 / 5 (f ) » (£) 1/5 p 1/5 

dominated (£) > (f) 1/2 p 1 /2 e -3/4 j > (»)V3 pl/ 3 £ -i/3 

regime (f ) » (f )- 1 p 1 /2 £ -3/4 (f j p i/4 £ -i/8 

(q/L) < 1 (a/L) < 1 

6 ~ qp 1 / 2 e" 1 '' 2 fo - ap 1/2 e-^ 2 

Panels in e 2 < p < e e 2 < p < e 

bending (f ) « (f ) 1/5 (£) > (f/V 8 

dominated (f) > (f ) 1/2 p 1/2 e" 3/4 (f) > (i^V^-i/a 

regime (f ) » (f ^V^- 9 / 4 (f ) » p 3 / 8 e " 3 / 8 

(q/L) < 1 (q/L) < 1 



TABLE II: The scaling properties of optimally efficient plates 
consisting of parallel spars supporting intervening panels. 
Spars and panels may be either in the bending or stretching 
dominated regimes. 





Spars in stretching 


Spars in bending 




dominated regime 


dominated regime 


Panels in 


p < e 5 / 2 


e 5/2 < p < e 2 


stretching 


Spars give 


(b) ~ Lp 3 / 4 e- 7 / 8 


dominated 


no gain 




regime 


in efficiency 


(») ~p-VV/ 8 


Panels in 


Not 


e 2 < p « e 


bending 


consistent 


(6) ~ Lp 5 / 8 e- 5 / 8 


dominated 


with 


(^) ^p3/8 £ -3/8 


regime 


conditions 


(«) ^pV8 e -V8 



V. FRACTAL DESIGN 
A. Two beams at right angles 



Because s <ti a, then over almost all their length, the 
spars must support a force per unit length given by / = 
2ap. 

We consider the four cases where the spars and panels 
are in the stretching and bending dominated regimes. 
There are six conditions which have to be applied in each 
case: 

(i) The panels are assumed to be at their most me- 
chanically efficient, which sets bo as a function of p and 
a. 

(ii) Whether the panels are in the bending or stretching 
dominated regime sets a constraint on the pressure. 

(iii) The spars are in the stretching dominated regime 
if ap 3> s 5 /L 4 and the bending dominated regime if ap <C 
s 5 /L 4 . 

(iv) The spars must not break (although we do not 
insist that they should be at this limit; they may be 
somewhat over-engineered in this sense) . 

(v) The spars must provide significant support, in the 
sense that the maximum out-of-plane displacement at the 
centre of each spar must be much less than the maximum 
out-of-plane displacement at the centre of each panel. 

(vi) There must be at least one spar in the system; in 
other words (a/L) < 1. 

These conditions are listed in Table |H 

It is a simple, but tedious matter to combine the con- 
ditions of Table [T] with Eq. to obtain the scaling 
properties of the optimally efficient solution in each case. 
The results are shown in Table HU 

Fig. 2] shows these behaviours on an o.o.m. plot, and 
(schematically) on the p — e plane. 



The parallel spars analysed above have provided con- 
siderable gains in efficiency over a uniform plate, but it 
is possible that we have not achieved the global optimal 
scaling. One possibility which remains is that we could 
allow spars to intersect. In this way, thicker spars may be 
used to partially support thinner spars, which together 
support panels. We can envisage ultimately a hierarchi- 
cal arrangement with spars of progressively thinner as- 
pect forming a fractal design, and supporting panels only 
at the smallest length scales. 

To proceed, we first investigate whether the principle of 
a thicker spar supporting a thinner spar can indeed lead 
to a gain in efficiency for a structure. Consider therefore 
a single spar which is in the bending dominated regime, 
carrying a (scaled) force / per unit length, and a point 
load F at the centre, which we represent by a scaled point 
force F = F/Y, The spar or beam must be considered 
in two sections, joined at the middle (x = 0), at which 
point the displacement w and the first two derivatives of 
w are continuous 

, . / w+(x) x G (0,L) , . 

W(X) — < _) N / r f\\ (25) 

w I w [x) x G (— L, 0) v ' 

Hence, if we have two beams at right angles with 
widths s\ and S2, joined at the centre, and each car- 
rying a force / per unit length, they will (in the absence 
of external point forces) apply equal and opposite forces 
F on each other, and their displacements will be given 
by the following expressions: 

«fM - j^j (*•-!»)■ ± jfg, if (.±f) 
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r1/2 - 




fore write the no-breakage condition on the first beam as 



log 10 £ 



FIG. 4: Top figure: order of magnitude behaviour for quan- 
tities as a function of pressure for an aperture spanned by 
a plate consisting of parallel spars supporting long panels. 
Bottom figure: schematic picture of the p — e plane, show- 
ing contours (solid lines) of constant number of spars in the 
structure. This number increases towards the bottom left of 
the figure. Region 'A' is where spars and panels are bend- 
ing dominated (e 2 <C p <C e). Region 'B' is where spars 
are bending dominated and panels are stretching dominated 
( e 5 / 2 p <g e 2 ^ an d region 'C is where there are no spars, 
and the now uniformly thick plate is stretching dominated 
(P«e 5/2 )- 



F = fL 



4 



(26) 



^|flgt^(*)| 



4/L 2 



G 2 s 3 (r 4 + 1) 



— max[|5r 4 -l|,(7r 4 + l)] <e 



> 



4/ 



L 3 ~ (r 4 + l)LG 2 e 



max[|5r 4 -l|,(7r 4 + l)] , (27) 



where r = s\/ S2- The result must be symmetric on swap- 
ping the two spars, so 



> 



if 



L 3 ~ (r- 4 + l)LG 2 e 



[|5r- 4 -l|,(7?- 



D] 



The total volume of material used to make these beams 
is therefore 



V = GL (sj + 4) 



(29) 



which we seek to minimise. 

Eqs. (|27|) and (|28|) have a trivial solution in which s\ = 
S2 and F = 0. The most efficient solution of this kind 
has both beams at the breaking limit, so 

8 X = s 2 = 2 4 / 3 L(//L) 1 / 3 e- 1 / 3 G- 2 / 3 , (30) 

and uses a volume of material given by 

V w 12.699L 3 (//L) 2 / 3 e- 2 / 3 G^ 1 / 3 . (31) 

However, this is not the global optimum, which in fact 
occurs when r ~ 1.22 and 



V w 12.652 J L 3 (//L) 2/3 e" 2/3 G" 1/3 , 



(32) 



and only the thinner of the two beams is at its breaking 
limit. 

We therefore see that although the gains are very mod- 
est, it is nevertheless possible to achieve greater efficiency 
through having a thicker beam or spar support a thinner. 
The next question to address is whether this principle can 
be extended to a hierarchical arrangement with thin spars 
supporting still thinner spars, until the very thinnest help 
to support panels to form a continuous plate with no 
holes. 



B. Hierarchical arrangement of spars 



where the point force F in Eq. (|26[) imposes the constraint 
that the beams must have the same displacement w(0) 
at their centre points. 

For this highly symmetrical case, the maxima in curva- 
ture occur either at the ends or the centres of the beams 
(for more general cases, there may be maxima between 
the positions of application of point loads). We can there- 



Consider the arrangement of spars as shown in Fig. [5] 
Fig. 03a) shows two spars at right angles, and we refer 
to this structure as "generation 1" . We choose the spars 
to have the same width si i (the first index referring to 
the "generation number", and the second to the type 
of spar present in this generation). Each spar carries a 
force per unit length /i = fo/2. In light of the previous 
calculation, we know that it is not optimal to have both 
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(a) 



(b)f 



TABLE III: Values for f3(n) and 7(71) (in Eq. (J33J)) which are 
found numerically to minimise V(n) (Eq. (|34|l from the text). 



n 


p{n) 


l{n) 


1 


1.260 




2 


1.126 


0.609 


3 


1.009 


0.644 


4 


1.078 


0.650 



(c) 



(d) 



the two parameters (3(n) and 7(71) denned by 



FIG. 5: Schematic of a fractal design for a plate spanning 
a square aperture, (a) is a "generation 1" structure, with 
two crossed spars supporting four panels. In (b), four thinner 
spars are added, which subdivides the remainder of the plate 
into sixteen panels. In (c) and (d) we show generations 3 and 
4. 



spars of the same width; but we ignore this fact in favour 
of simplicity of analysis. 

Fig [5fb) shows a "generation 2" structure, where we 
have two thick spars, just as for generation 1, but their 
widths are now S2,i which may be different from s^i. In 
addition, there are four thinner spars, each with width 
S2,2- All the spars now carry a force per unit length of 
h = /o/4 - chosen to approximate the loading produced 
by the smaller panels between the spars which will even- 
tually be present in the structure. Fig.(5[c) and (d) show 
generations 3 and 4 respectively, and we envisage going 
to large generation numbers n for certain regions of the 
p — e plane. 

For generation n, we will have two spars of width s n ,i, 
four thinner spars of width s n> 2 and so on, culminat- 
ing with 2™ of the thinnest spars, which each have a 
width s n _ n . The total number of spars in the genera- 
tion n structure is therefore 2(2™ — 1), and each carries 
a force per unit length (in addition to the point forces) 
of /„ = 2~"/o. Again, this is chosen to mimic the load- 
ing on the spars produced by supported panels, so that 
fo w pL/A (although the loading in this final structure 
will not be exactly uniform). 

We now analyse the first few generations numerically, 
ignoring the possible torsional loadings that spars can 
exert upon one another in asymmetrical configurations. 
We expect this complication to have minimal effects when 
G < 1 , and so we now specialise to the case G = 1 . 

For generation n, we choose a structure described by 



L 



1/3 



,-1/3 



f3(n) [ 7 (n)] 



q-l 



(33) 



We choose (3(n) and 7(71) so that all the spars satisfy 
the no-breaking condition, and the amount of material V 
that is used, is minimised. Adding up the total volume 
of material in the spars gives 



V(n) = 2L 3 



2/3 



c-2/3fl2 



[27 2 (n)]' 



1 



27 2 (n) - 1 



(34) 

Based on experience from the previous section, we ex- 
pect that the most efficient structure of this kind will 
satisfy the breakage condition as an equality only for the 
thinnest spars. We also expect 7(n) > 0.5 (the thicker 
spars will not act as completely rigid supports for the 
thinner spars), but we hope that a hierarchical structure 
is still achieved, which corresponds to "f(n) < 1. 

Table [TTT1 shows the first few values for /3(n) and ~f(n) 
which arc found from numerical simulations to minimise 
V{n). The tentative conclusion is that both quantities 
tend to fixed, order 1 values as n increases. 

If now we assume that P(n) and "f{n) do indeed tend 
to fixed limits (3 and 7 as n — > 00, then provided the 
spars remain much narrower than the panels, and the 
spars remain in the bending dominated regime, we can 
calculate the average thickness of the plate. It will be a 
sum of two terms; one from the hierarchical structure of 
spars, and one from the panels between (which may be 
in either the stretching or bending dominated regimes). 
For an order of magnitude calculation, the pre-factors are 
unimportant and we expect square panels to exhibit the 
same scaling as the long thin panels analysed above. By 
using Eqs. (j3~4"|) and (fT5)) we obtain in the large n limit, 
and assuming (3 is an order 1 quantity 

<6) / (2 7 2 ) V /3 e~ 2/3 + 2->e- 3 / 2 p«e 2 

L ~ I (2 7 2 )> 2/ V 2/3 + 2- V /2 e~ 1/2 e 2 «P«e ' 

If we minimise this as a function of n, we find 



2 ^ \ In 2 In 2 



p 121n(2 7 ) e 121n(2 7 ) £ '«j,<^ ( 



(A 

L 



pi 1 6 ln(2~f) )g( 2 + 12 J ln*(2 7 ) ) 
5 ( T + 12 inC2 7 ) ) g ( - 2 - 12 !S(S 7 ) ) 



p < e 

e 2 < p < e 



(36) 



. (37) 
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where for each of Eqs. (|36|) and (|37|) . the upper expression 
corresponds to panels in the stretching regime and the 
lower to panels in the bending regime. 

We therefore find that the hierarchical arrangement 
is more efficient than a uniform plate, provided that 
£ 5/2 <g- p ^ £j but for smaller values of p, a plate of 
uniform thickness is most efficient. This is precisely the 
range over which parallel spars are more efficient than 
uniformly thick plates, so the next question is whether a 
fractal structure is more efficient than the parallel spars? 

By comparing Eq. (j3"Tj) with the two expressions for (b) 
in Table |TT] we find that there is a critical value 7 C of 7, 
above which the parallel arrangement is more efficient, 
and below which fractal structures are preferred. The 
critical value is 

7c = 2^/ 3 w 0.7937. (38) 

Based on the numerical calculations of Table 11111 it 
appears that 7 ss 0.65 < 7c , and so (provided e <C 1 in the 
mathematical sense), then the hierarchical structure is 
always more efficient over the entire range e 5 / 2 «jj«6. 

We note however that for small, but finite values of e, 
the pre- factors are no longer negligible, and it remains a 
possibility that there is a region of the p — e plane not too 
far from the origin (log e, log p) = (0,0), where a parallel 
arrangement of spars is preferred. 

Lastly, we note that for the fractal structure, the width 
of the thinnest spars is given by 

s n ,n J p(^ + 6i 1 °(?T) ) e (~3- mzfcj ) p<e 2 

L \ p(^~ 121n(2 7 ) ) e (~l+ 121n(2 T ) ) £ 2 < J) < £ ' 

(39) 

so s n , n < « = 2~ n L for e D / 2 < p < e and 7 G (0.5, 1), 
which justifies Eq. (|35| for (b). 



VI. CONCLUSIONS 

We have constructed a simple problem in elasticity the- 
ory, in which it is possible to frame the question of opti- 
mal mechanical efficiency in an easily approachable man- 
ner. We find that it is necessary to consider all regions of 
the p— e plane to tackle this problem (where p represents 
the applied pressure and e the material's brittleness) and 
have used order-of-magnitude plots as a concise way to 
present the assymptotic behaviour under the assumption 
e<l. 

After considering various possible forms for the plate 
in the problem, we discover regions of the p — e plane 
where a fractal design is the most mechanically efficient 
structure we have been able to find. This is a structure 
where thicker spars act as partial support for thinner 
spars, and so on until the thinnest spars (together with 
all the others) support panels, which form a continuous 
plate. 

Our analysis has focused on the order of magnitudes of 
quantities; this greatly simplifies the algebra, but leaves 
open questions about the behaviour for small, but finite 
values of p and e. 
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